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SUMMARY 


Solutions  to  the  governing  equations  of  a  spur  gear  transmission  model, 
developed  in  NASA  TM-100180  (AVSCOM  TM-87-C-2)  are  presented.  Factors  affect¬ 
ing  the  dynamic  load  are  identified.  It  is  found  that  the  dynamic  load 
increases  with  operating  speed  up  to  a  system  natural  frequency.  At  operating 
speeds  beyond  the  natural  frequency  the  dynamic  load  decreases  dramatically. 
Also,  it  is  found  that  the  transmitted  load  and  shaft  inertia  have  little 
effect  upon  the  total  dynamic  load.  Damping  and  friction  decrease  the  dynamic 
load.  Finally,  tooth  stiffness  has  a  significant  effect  upon  dynamic  loading: 
the  higher  the  stiffness,  the  lower  the  dynamic  loading.  Also,  the  higher  the 
stiffness  the  higher  the  rotating  speed  required  for  peak  dynamic  response. 


INTRODUCTION 


The  development  of  a  simple  parallel  shaft  spur  gear  transmission  model 
with  its  dynamic  differential  equations  and  solution  procedures  were  presented 
in  reference  1.  Various  parameters  such  as  inertia,  stiffness,  friction  and 
damping  were  included  in  the  governing  equations  for  further  study  (ref.  1). 


The  purpose  of  this  report  is  to  determine  the  effect  of  these  parameters  i>  For 
on  gear  dynamic  load.  The  dynamic  load  of  a  gear  transmission  can  be  found  by  ,vI 
solving  the  governing  equation.  The  solution  is  known  as  the  dynamic  motion 
of  the  gear  transmission.  This  dynamic  motion  can  then  be  substituted  into 
other  analytical  formula  and  solved  for  gear  dynamic  loads.  ... 
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A  model  of  the  transmission  is  depicted  in  figure  1.  The  governing  equations 
are: 

JM0M  +  Csl<0M  ~  01}  +  Ksl(6M  ~  Gl)  =  TM  ( 

J101  +  Csl(0r  0M}  +  Ksl(9l  “  0M>  +  Cg(t)  [Rbl0  ”  Rb202]  (, 

+  Kg(t)  CRbl(Rbler  Rb2e2)]  -  Tfl(t) 

J202  +  Cs2(02  ~  01}  +  Ks2 (02  "  01 }  +  Cg(t)  CRb202  “  Rbl0l ]  (. 

+  Kg<t)  CRb2<Rb202  -  Rblo,)]  =  rf2(t) 


JL0L  +  Cs2  <0L  ”  °2)  +  Ks2<0L  "  °2)  =  ”TL  (4) 

where  J^,  Jj,  J2 .  and  J(_  represent  the  mass  moments  of  inertia  of  the  motor, 
the  gears,  and  the  load;  Csi,  Cs2,  and  Cg(t)  are  damping  coefficients  of  the 
shafts  and  the  gears;  Ks i ,  Ks2,  and  Kg(t)  are  stiffness  of  the  shafts  and  the 
gears;  T^,  T|_,  Tf|(t),  and  Tf2(t)  are  motor  and  load  torques  and  frictional 
torques  on  the  gears;  Rb]  and  Rb2  are  base  circle  radii  of  the  gears;  t  is 
time,  and  the  overdots  indicate  time  differentiation. 

In  this  report  we  present  the  results  of  numerical  solutions  of  equa¬ 
tions  (1)  to  (4)  for  a  typical  transmission  system.  A  flow  chart  outlining 
the  numerical  procedure  is  presented.  Natural  frequencies  are  determined. 

The  dynamic  load  is  determined.  Finally,  the  dynamic  factor  defined  as  the 
ratio  of  the  dynamic  load  to  the  static  load,  is  determined.'  The  results  are 
calculated  as  functions  of  the  rotating  speed  and  roll  angle  for  a  variety  of 
damping  and  stiffness  conditions. 


NOMENCLATURE 

damping  coefficient  of  gear  tooth  mesh,  N-sec  (lb-sec) 
damping  coefficient  of  shaft,  N-m-sec  (in. -lb-sec) 
polar  moment  of  inertia  of  load,  kg-m2  (in. -lb-sec2) 
polar  moment  of  inertia  of  motor,  kg-m2  (in. -lb-sec2) 
polar  moment  of  inertia  of  gear  1,  kg-m2  (in. -lb-sec2) 


'The  term  "dynamic  factor"  or  "dynamic  load  factor"  has  been  used  inconsis¬ 
tently  in  the  literature.  The  American  Gear  Manufacturer 1 s  Association 
(AGMA)  dynamic  factor  Kv  is  used  as  a  strength  reduction  factor  and  is 
defined  as  the  maximum  static  load  divided  by  the  maximum  dynamic  load.  This 
paper  will  follow  the  ISO  convention  which  uses  the  dynamic  factor  Kb  as  a 
load/stress  increasing  factor.  Therefore,  Kb  =  1/KV. 
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polar  moment  of  inertia  of  gear  2,  kg-m2  (in. lb-sec2) 


dynamic  factor 


stiffness  of  gear  tooth,  N-m/rad  (in. -lb/rad) 


stiffness  of  shaft,  N-m/rad  (in. -lb/rad) 


AGMA  dynamic  factor,  Kv  =  l/K^ 


base  radius,  mm  (in.) 


pitch  radius,  mm  (in. ) 


torque  on  load,  N-m  (in. -lb) 


torque  on  motor,  N-m  (in. -lb) 


torque  on  gear  1,  N-m  (in. -lb) 


torque  on  gear  2,  N-m  (in. -lb) 


applied  load,  N  (lb) 


angular  displacement,  rad 


angular  velocity,  rad/sec 


angular  acceleration,  rad/sec2 


natural  frequency,  Hz 


damping  ratio 


PROCEDURE 


Figure  2  presents  a  flow  chart  of  the  computational  procedure  used  in  the 
parameter  study.  The  procedure  is  the  same  as  that  outlined  at  the  end  of 
reference  1 . 


In  conducting  the  analysis  it  is  useful  to  compare  the  locations  of  the 
peak  dynamic  loads  with  the  locations  of  the  system  natural  frequencies  (or 
critical  speeds).  The  natural  frequencies  themselves  may  be  obtained  by  exam¬ 
ining  the  undamped  equations  of  motion.  These  equations  may  be  written  in  the 
matrix  form: 


CJ]  [9]  +  CK]  [0]  =  [0] 


where  the  inertia  matrix  [J]  is 


JM  0  0  0 


[J]  = 


0^00 


0  0  j2  0 


0  0  0  J, 


‘.'w  vv -r  ••  /  V  V.  Sm  «\  n*  s.  .  •r. 


and  the 

stiffness 

matrix 

[K]  i  s 

'  Ksl 

'Si 

0 

0 

[K] 

. 

-Ksl 

Ksl  + 

Vav/b, 

-(K  )  R..R... 

g  avg  bl  b2 

0 

0 

_(Kg 

)  R,  ,  R,  „ 
avg  bl  b2 

Si  ♦  (Kg>avgR2b2 

-Ks2 

0 

0 

Ks2 

where  (kg)aVg  represents  the  average  value  of  the  gear  mesh  stiffness.  It 
is  taken  as  tne  sum  of  the  discrete  tooth  stiffness  values  of  a  mesh  cycle 
divided  by  the  number  of  mesh  positions  in  the  cycle. 

In  the  parameter  study  the  system  had  identical  gears  with  the  following 
properties: 


Number  of  teeth  . 

Module,  mm  . 

Pi tch  diameter ,  m  ( i n . )  . 

Pressure  angle,  deg  . 

Appl i ed  load ,  N  ( 1 b)  . 

Face  width ,  m  ( i n . )  . 

Moment  of  inertia,  kg  m2  (in. -lb  sec2)  .  . 
Average  tooth  stiffness,  N  m/rad  (lb-in. /rad) 
Damping  ratio  . 


.  36 

3.18  (8  di ametra 1  pitch) 

.  0.1143  (4.5) 

.  20 

.  2670  (600) 

.  0.0254  (1 .0) 

.  3.3323xlO-3  (0.02947) 

3.991xl05  (3.5355xl06) 
.  0.10 


The  shaft  stiffness  and  inertias  were: 


Shaft  stiffness,  N  m/rad  (in. -lb/rad)  .  1138.17  (10081) 

Load  inertia  and  motor  inertia, 

kg  m2  (in. -lb  sec2)  .  9 . 989x 1 0-3  (each)  (0.08841) 


RESULTS 


Using  the  aforementioned  data  in  the  gear  system  model  shown  in  figure  1, 
the  natural  frequencies  of  the  four  degrees  of  freedom  model  were  found  to  be 


wnl  =  0  (rigid  body  mode) 

wn2  =  1.49  Hz 

or  .  =  2.99  Hz 
n3 

un4  =  Hz 

(8) 


The  first  three  natural  frequencies  are  well  below  tooth  meshing  frequencies 
and  are  therefore  not  of  interest  in  this  analysis  (although  they  can  still  be 
excited  by  shaft  eccentricity  which  is  not  modeled  here).  The  fourth  natural 
frequency  matches  tooth  meshing  frequency  at  the  critical  speed  of  8688  rpm 
which  is  within  the  operating  range  of  the  gears. 


Figure  3  shows  the  variation  of  dynamic  load  response  for  a  pair  of  teeth 
as  a  function  of  roll  angle.  At  speeds  much  lower  than  the  critical  speed, 
the  dynamic  load  response  is  basically  a  static  load  sharing  in  phase  with  the 
stiffness  change,  superimposed  with  an  oscillatory  load  at  a  frequency  corre¬ 
sponding  to  the  natural  frequency. 
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At  higher  speed,  close  to  the  critical  speed,  the  dynamic  load  variation 
becomes  so  abrupt  that  it  produces  tooth  separation.  The  peak  dynamic  load  is 
much  higher  than  the  static  load  and  is  very  likely  to  be  a  source  of  gear 
noise  and  early  surface  fatigue. 

Figure  4  shows  the  dynamic  factor  iCq  as  a  function  of  operating  speed. 
Prominent  peaks  (resonances)  may  be  seen  at  speeds  of  7650  and  4200  rpm.  The 
larger  peak  7650  rpm  occurs  at  38  percent  of  the  calculated  critical  speed. 

The  experimenta1  work  by  Kubo  (ref.  6)  reported  a  similar  result  that  the  crit 
ical  speed  was  found  at  about  90  percent  of  the  calculated  critical  speed. 

The  second  peak  at  4200  rpm  is  a  nonlinear  effect  of  the  time  varying  tooth 
stiffness  known  as  the  parametric  resonance.  This  parametric  resonance  fre¬ 
quency  occurs  at  about  one-half  of  the  critical  speed  (ref.  2).  For  speeds 
above  the  critical  speed,  the  dynamic  response  decreases  steadily  in  the  same 
manner  as  with  elementary  vibrating  systems. 

Figure  5  shows  a  three-dimensional  representation  of  the  system  dynamic 
response.  The  horizontal  axis  represents  the  operating  speed,  and  the  contact 
position  along  the  tooth  profile.  The  total  number  of  contact  positions  is 
121.  The  vertical  axis  is  the  dynamic  factor. 

Figure  6  presents  a  contour  plot  of  the  system  dynamic  response.  The 
shaded  areas  represent  regions  where  tooth  separation  occurs.  They  are  locat¬ 
ed  in  the  double  contact  regions.  At  near  resonance  speeds  the  vibration 
amplitude  exceeds  the  deflection  of  the  meshing  teeth,  thus  inducing  tooth 
separation . 

As  the  speed  increases,  the  dynamic  response  also  shows  a  phase  shift 
towards  the  higher  numbered  contact  positions.  This  phenomenon  can  be  seen  by 
noting  that  at  speeds  from  600  to  4200  rpm  (one-half  subharmonic),  the  maximum 
dynamic  load  occurs  at  double-  to  single-contact  transitions  (position  51)  and 
gradually  changes  to  single-  to  double-contact  transitions  (position  75).  At 
speeds  between  one-half  subharmonic  and  resonance  the  maximum  peak  stays  near 
the  single-  to  double-contact  transitions.  After  the  speed  passes  resonance, 
the  major  dynamic  peak  moves  again  towards  higher  contact  positions  on  tooth 
profile  with  increasing  speed. 

Figure  7  shows  the  dynamic  factor  as  a  function  of  the  transmitted  loads 
for  three  different  speeds.  There  is  only  a  small  decrease  in  dynamic  factor 
with  increased  applied  load. 

Figure  8  shows  the  effect  of  damping  on  the  dynamic  load.  It  is  seen 
that  damping  has  its  greatest  effect  near  resonance  frequencies. 

Changes  in  shaft  stiffness  have  a  minor  effect  on  the  system  dynamic 
response.  However  changes  of  tooth  stiffness  have  a  major  effect  on  the 
response.  Figure  9  shows  that  the  higher  the  tooth  stiffness  the  lower  the 
dynamic  response  (dynamic  factor).  This  is  consistent  with  observations  that 
as  the  tooth  stiffness  increases  the  effect  of  gear  mass  on  the  system  dynam¬ 
ics  is  reduced.  Figure  9  also  shows  that  system  resonance  frequencies  are 
increased  as  the  tooth  stiffness  increases.  This  is  a  potentially  useful 
effect  for  the  design  of  gear  systems. 

The  effect  of  shaft  mass  is  assumed  small  compared  to  that  of  the  gears. 
Figure  10  shows  that  as  the  gear  inertia  is  reduced,  the  dynamic  response  is 
al so  reduced . 


For  gears  with  different  diametral  pitches,  the  dynamic  response  is  dif¬ 
ferent  due  to  the  change  in  contact  ratio.  Gears  with  a  finer  pitch  have  a 
higher  contact  ratio.  Since  the  contact  ratio  is  a  measure  of  the  duration  of 
the  load  being  shared  by  more  than  one  pair  of  teeth,  it  has  a  significant 
effect  on  the  system  dynamic  response. 

Figure  11  shows  a  comparison  between  gears  having  different  diametral 
pitches.  The  finer  pitch  gears,  having  a  higher  contact  ratio,  have  a  smaller 
dynamic  loaa  than  the  coarser  pitch  gears. 


DISCUSSION 


In  1927,  A. A.  Ross  (ref.  3)  introduced  the  following  empirical  formula 
for  the  dynamic  factor  kv: 


k 


v 


78 

(78  +  -y/v) 


(9) 


where  v  is  the  pitch  line  speed  measured  in  ft/min.  This  expression  received 
acceptance  as  a  standard  factor  used  by  the  American  Gear  Manufacturer's 
Association  (AGMA).  In  1959,  a  similar  factor  for  use  with  higher  precision 
gears  was  introduced  by  Wellauer  (ref.  4) 


k 


v 


78  +  «y/v 


(  10) 


Equations  (9)  and  (10)  are  recognized  as  being  conservative  when  applied 
with  very  high  precision  gears.  They  are  thought  to  predict  dynamic  loads 
which  are  higher  than  the  actual  loads. 

Buckingham  (ref.  5)  has  also  developed  an  expression  for  the  dynamic  load 
in  terms  of  the  pitch  line  speed  and  the  applied  load.  His  formula  is 

W  ,  =  W  +  [f  (2fh-  f  >]'/l  (11) 

d  aba 


where  is  the  dynamic  load,  W  is  the  applied  load  and  the  factors  fa 

and  fp  are 

f  =  f.f  /( f,  +  f  )  and  f,  =  0.0000555  EF  ♦  W  (12) 

a  b  c  b  c  b 

where  F  is  the  face  width,  E  is  the  elastic  constant,  and  fc  is 

f  =  0.00120  C ( R 1  ♦  R ^ ) / R | R 2 Imv^  (13) 

where  Ri  and  R?  are  the  pitch  radii  of  the  gears  and  m  is  their  effec¬ 
tive  mass.  (In  these  expressions  the  units  are  in  pounds  and  inches  except 
for  the  pitch  line  speed  which  is  measured  in  ft/min.) 


5 


I.TV 


Kubo  (ref.  6)  measured  static  and  dynamic  gear  stresses  on  several 
high-preci sion  spur  gear  systems.  Kubo  expressed  the  dynamic  factor  as  the 
ratio  of  maximum  dynamic  to  the  maximum  static  stress.  Since  stress  is  propor¬ 
tional  to  load,  Kubo's  definition  of  dynamic  factor  is  identical  to  that  used 
here.  Figure  12  shows  a  comparison  of  (1/KV)  from  the  AGMA  high-precision  for¬ 
mula  (eq.  (10)),  (W(j/W)  from  Buckingham's  formula,  Kubo's  results,  and  the 
results  of  the  computer  simulation  for  an  identical  spur  gear  pair  with  4  mm 
module,  25  teeth,  20°  pressure  angle,  15  mm  face  width,  131.5  kN/m  applied 
load,  and  207  GPa  Young's  modulus. 

There  is  a  good  agreement  between  Kubo's  result  and  the  computer 
simulation. 


CONCLUSIONS 

From  the  foregoing  results  several  conclusions  may  be  stated: 

1.  For  accurately  made  gears,  the  dynamic  load  is  significantly  affected 
by  the  contact  ratio:  for  increased  diametral  pitch,  that  is,  for  high  contact 
ratio  gears,  the  dynamic  load  is  decreased. 

2.  The  tooth  stiffness  has  a  significant  effect  upon  the  dynamic  load: 
the  higher  the  stiffness  the  lower  the  dynamic  load.  Also,  the  higher  the 
stiffness  the  higher  the  critical  speed  for  peak  response. 

3.  The  dynamic  load  generally  increases  with  the  operating  speed  until  a 
resonance  is  reached.  The  dynamic  load  decreases  rapidly  beyond  the  resonance. 

4.  Damping  and  friction  decrease  the  dynamic  load  with  the  most  dramatic 
effects  occurring  near  the  resonance  frequencies. 

5.  The  applied  load  has  a  relatively  minor  effect  upon  the  dynamic  factor. 

6.  Tooth  separation  —  leading  to  impact  --  occurs  in  the  double  tooth 
contact  region  since  the  deflections  are  smallest  in  that  region. 

7.  The  dynamic  factor  is  largest  for  contact  points  near  the  tooth  tip. 

8.  The  dynamic  factor  decreases  with  decreases  in  the  gear  body  inertia. 
Shaft  moment  of  inertia  has  a  minima!  effect  upon  the  dynamic  factor. 
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16  Abstract 


Solutions  to  the  governing  equations  of  a  spur  gear  transmission  model,  developed  in  NASA  TM- 100180 
(AVSCOM  TM-87-C-2),  are  presented.  Factors  affecting  the  dynamic  load  are  identified.  It  is  found  that  the 
dynamic  load  increases  with  operating  speed  up  to  a  system  natural  frequency.  At  operating  speeds  beyond  the 
natural  frequency  the  dynamic  load  decreases  dramatically.  Also,  it  is  found  that  the  applied  load  and  shaft  inertia 
have  little  effect  upon  the  dynamic  load.  Damping  and  friction  decrease  the  dynamic  load.  Finally,  tooth  stiffness 
has  a  significant  effect  upon  dynamic  loading:  the  higher  the  stiffness,  the  lower  the  dynamic  loading.  Also,  the 
higher  the  stiffness  the  higher  the  rotating  speed  required  for  dynamic  reponse. 
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